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Abstract 

We compute the electric screening mass in lattice QCD with Wilson fermions 
at finite temperature and chemical potential to one-loop order, and show that 
lattice artefacts arising from a finite lattice spacing result in an enhancement of 
the screening mass as compared to the continuum. We discuss the magnitude 
of this enhancement as a function of the temperature and chemical potential 
for lattices with different number of lattice sites in the temporal direction that 
can be implemented in lattice simulations. Most of the enhancement is found 
to be due to the fermion loop contribution. 

PACS numbers: 11.15.H, 12.38. G 
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1 Introduction 



An important feature of finite temperature QCD is the generation of an electric and 
magnetic screening mass which play an important role in controlling the infrared 
behaviour of the theory. The electric screening mass leads to a Debye screened static 
quark- antiquark potential and is given for SU(N) with Nf quark flavours, and for 
vanishing chemical potential and quark mass, in leading order perturbation theory 
by '^ei — ^{N + Nf /2)T'^ Renormalization group improved perturbation theory 
tells us that the effective coupling is a function of the temperature, and decreases 
with increasing temperature. This suggests that at sufficiently high temperatures, 
above the deconfining phase transition, the screening mass may be computed in 
perturbation theory. Because of the singular infrared behaviour of the perturbative 
series, however, the computation of the next to leading order contribution requires 
a resummation of infrared divergent diagrams which turns out to be sensitive to 
the magnetic screening mass 0. This mass vanishes in lowest order perturbation 
theory and is expected to be of 0{g'^T). The coefficient multiplying g'^T turns out 
however to be incalculable 0. Making use of an improved perturbation theory 
proposed by Braaten and Pisarski [Q, which resums hard thermal loops, and of a 
gauge invariant definition of the electric screening mass Rebhan has calculated 
the 0{g^T'^) corrections to the non-abelian screening mass-squared and has shown 
that next to leading order contributions give rise to an enhancement p|. 

The lattice formulation of QCD allows one to determine the electric screen- 
ing mass non-perturbatively. The screening mass is extracted from correlators of 
Polyakov loops [6-9], or from the long distance behaviour of the gluon propaga- 
tor [jlO|. For small quark-antiquark separations lattice perturbation theory for the 
Polyakov loop correlation function is expected to describe the Monte Carlo data, 
since for a finite lattice volume one is not confronted with the infrared problems 
encountered in thermal perturbation theory. This has been checked in for the 
SU(3) gauge theory by taking careful account of finite size effects and, in particular, 
of the zero momentum modes which do not allow one to take the thermodynami- 
cal limit for fixed coupling, as one would do in standard perturbation theory. For 
larger quark-antiquark separations, beyond the "perturbative horizon", the colour 
averaged potential is expected to have a Debye screened form. Monte Carlo simula- 
tions confirm this screening picture [6-9]. In the case of pure SU(2) and SU(3) gauge 
theories the electric screening mass, when determined from Polyakov loop correlation 
functions, is found to be about 10% larger [|, 0] than the leading order perturbative 
result if the temperature dependent coupling constant is determined from Polyakov 
loop correlators in the perturbative region. As was pointed out by Rebhan, [Q , such 
an enhancement could also be expected if next to leading order corrections to the 
continuum screening mass are taken into account through resummed perturbation 
theory. A quantitative comparison with the results obtained in the above simulations 
is however very difficult and has, to our knowledge, not been carried out so far. In 



contrast to the work of ref. [o, 0, the electric screening mass as extracted in ref. [10 



2 



from the gluon propagator in the Landau gauge was found to deviate strongly form 
the leading order perturbative result. 

In comparing the Monte Carlo data for the electric screening mass with leading 
order, or resummed, perturbation theory it is important to have an estimate of the 
size of lattice artefacts to be expected from a finite lattice spacing. To obtain such an 
estimate we compute the electric screening mass in one-loop order on the lattice, at 
finite temperature and chemical potential, in the infinite volume limit, and compare 
it with the continuum. For the case of QED with naive fermions the screening 



mass has been calculated by Pietig |]12|. The screening mass is calculated from 



the zero-momentum limit of the 44-component of the vacuum polarization tensor 
evaluated for vanishing Matsubara frequency. In one-loop order this definition of 
the screening mass is gauge invariant, and consistent with the more general gauge 
invariant definition given in ref . , where the screening mass is determined from the 
position of the pole of the gluon propagator for vanishing Matsubara frequency. 

The paper is organized as follows: In the following section we calculate the electric 
screening mass for QCD in one-loop order for the case of Wilson fermions. The 
Feynman rules and frequency summation formulae required for the computation are 
relegated to two Appendices. As we shall see, the resulting integral expression has a 
very transparent form. In section 3 we then evaluate the momentum integrals for the 
screening mass numerically and compare the results with the continuum. We show 
that the lattice artefacts due to a finite lattice spacing give rise to an enhancement of 
the screening mass as compared to the continuum. We discuss the magnitude of this 
enhancement as a function of the temperature and chemical potential for lattices with 
different number of lattice sites in the temporal direction which can be implemented 
in numerical simulations. Most of the enhancement is found to be due to the fermion 
loop contribution. Section 4 contains a summary of our results. 



2 Electric Screening Mass in One Loop Order 

In this section we compute the electric screening mass in lattice QCD to one loop 
order from the zero momentum limit of the 44-component of the vacuum polariza- 
tion tensor evaluated for vanishing Matsubara frequency. The Feynman diagrams 
contributing in this order are shown in fig. 1. While diagrams (a,c,d,e) have a con- 
tinuum analog, the remaining diagrams, required by gauge invariance on the lattice, 
do not possess a continuum counterpart. The finite temperature, finite chemical po- 
tential lattice Feynman rules are obtained from the T = fi = rules by replacing the 
fourth component of the fermion and boson momenta by a)^ + ^''^d a)/, respec- 
tively, where u)/ = {lj^ = -^^^i^), with £ G Z, are the Matsubara frequencies 

for bosons (fermions), and f3 is the inverse temperature. Quantities with a "hat" are 
always understood to be measured in lattice units. Furthermore, integrals over the 
fourth component of momenta at zero temperature are replaced at finite temperature 
by sums over Matsubara frequencies in the interval [— |, | — 1], where we have taken 
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(a) (b) (c) (d) 



(e) (f) (g) 

Figure 1: Feynman diagrams contributing to the vacuum polarization in one loop 
order on the lattice. 

P to be even. The Feynman rules are collected in Appendix A. The relevant formulae 
for carrying out the sums over Matsubara frequencies are derived in Appendix B. 

The contributions of diagrams (a-g) to the vacuum polarization tensor are diag- 
onal in colour space, 

ri(f/)^^(a;/, k) = SABflfj^\u;i, I) . (2.1) 
The electric screening mass (in lattice units) is then defined by 

m^^^hmng-'^Ho,^) . (2.2) 

— * 

In the following we first consider the contributions to 1144''^^ (0,^) coming from the 
fermion loops, i.e., diagrams (a) and (b). 

i) Contribution of diagram (a) 

A straight forward application of the finite temperature, finite chemical potential 
lattice Feynman rules yields 

nr/^(0, "k\a) = -^A E r ^J^^V^'^'-'^l k) , (2.3a) 



where 



with 



f) . (2.3b) 



11 . ^ 

+ . ^ r > (2-3c) 



[1 + M(p)] + 
4 



Q 

M{p) 



[l + Mip)][l + MiP + k)] 

1 + ^ sin pj sin(p + k)j , 
j 

m + 2 ^ sin^ ^ , 



(2.3d) 
(2.3e) 

(2.3f) 



and Nf is the number of quark-flavours. The position of the poles of f^°'\z;p, k) are 
given by 



e^; Z2 = e 



with 



= ^^(P), 

= £{p + A;) 



(2.4a) 



where 



and 



Note that 



8{q) = In 
Kiq) = 1 



K{q) + ^K\q) - 1 



arcosh_R'(g) 



2[1 + M(g)] 



E(g)= /5:sin2g, +M2(g) . 



< > yj 



(2.4b) 
(2.4c) 

(2.4d) 
(2.5) 



while rj, C, and Q are invariant under the transformation p — ^ —p — k. This will be 
important further below. The frequency sum can be performed making use of eq. 
4]) derived in Appendix B. One finds that 



1 



1 



e/3(</.+A) _^ 1 ^-i3{4,-ii) I 
1 1 



where 



cosh 2(j) — T] cosh + 
sinh 0(cosh — cosh ip) 



(2.6a) 
(2.6b) 
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To obtain 1144'^-' (0; k)(^a), we must integrate this expression over p, with pj G [— vr, 

Noting that r], ^ and Q are invariant under the transformation p —p — k, and 
making use of ( p.5| ), as well as of the fact that the integrand in ( |2.3a| ) is a periodic 
function in pj and kj, we can combine the last two contributions on the r.h.s. of 
( p.6a| ) and obtain 

-Nfg' £^-0^hi(l),^lj,7],^,g)[ripni(i))+fiFDi<l>)] ■ (2.Ta) 

1 - 1 

^FD{(f)) = A., , ^ , Vpj^i^) = , (2.7b) 

are the lattice Fermi-Dirac distribution functions for particles and antiparticles. 

a) Contribution of diagram (b) 

We next compute the contribution to 1144'^'' (0, A;) of the Feynman diagram (b) 
depicted in Fig. 1. This diagram has no analog in the continuum and is given by 





IT 



where 



ni^'^^O, ^)(,) = -Nfg^- E r ^f^ie^^^^'^^J) (2.8a) 

P p-i •'^^ ^'^'^> 
*~ 2 



where 

f'\-3) = - "'"'f + \ , (2.8b) 
{z - zi)[z - Z2) 

p = , (2.8c) 

1 + M{p) 

and where zi and Z2 have been defined in (2.4a-d). Making again use of the frequency 
summation formula (IB. 41), one verifies that 



nir(o.f).,^-A',//;|fj(coth,--^-i) 

Combining this expression with (|2.7a] ) one finds that 

nl^'^)(oi) = nir^(o,?) + Ar^/£^i/(0,^,p,^ , 

(2.10a) 
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where 



H{(j), ^, p, r/, e, g) = coth - - h{<p, ^, 7], e, g) (2.10b) 



and 

nlr^(0,^) = -iV;/£^if(0,^,P,r/,e,^) (2.10c) 

is the T = /i = contribution. As we now show, n44"^''(0, A;) vanishes in the hmit 
/c — i> 0, and hence does not contribute to the screening mass. 



Consider the function h{(j),ilj,ri,^,g) defined in (|2.6b|) . It is singular for — 0, 
since in this hmit ip ^ (p. The singularity is however integrable.Q This can be seen 
as follows. Since according to (|2.5|), and the statement following it 



/i(0, i:, 7^, e, g) hiij, 0, r/, e, g) (2.11) 



we can also write ( [^.lOc] ) in the form 



fiir)(o,^) = -N f 9' f_^^A coth 



(2n)^ \ sinh( 
(2^ 



+ iNfg' 7|^^('^' ^' ^' ^' ^) ' (2.12a) 



where 



/i(0, r/, e, 6^) = /i(0, V-, ^, e, 6^) + Kij, 0, r/, 6^) . (2.12b) 



Although each term on the rhs of ( p.l2b| ) is singular for /c — ^ {ip ^ (f)), the sum 



possesses a finite limit. Thus setting tp = (p + e and taking the limit A; ^ (e ^ 0), 
one verifies that 

lim/i(0,V^,r7,e,e^) = 2 (coth0- -4- ) • (2.13) 
k^o V smh0; 



From ( |2.12a| ) we therefore conclude that 



limflir^(0,A;) = . 

This result is not unexpected, since for vanishing temperature and chemical potential 
it is well known in the continuum formulation, that Lorentz and gauge invariance 
protects the gluon from acquiring a mass. The screening mass is therefore determined 
by the finite temperature (f.T.), finite chemical potential contribution, given by the 

integral in ( |2.10a| ). By making again use of the fact that (p *->■ tp, when p —p — k, 

^ We follow here and in the following a technique used in ref. , where the author has calculated 
the screening mass for naive fermions in lattice QED to one loop order 



7 



while rj, ^ and Q remain invariant under this change of variables, we can write this 
contribution in the form 

+ h{^,<P,r],^,g)[fiFDW +flFDm} ■ (2.14) 
Consider the second integral. It can be rewritten as follows 

where h has been defined in ( |2.12b| ), and 

Mfd{<P, ^) = [f]FD{4>) - flFoW] + [flFDi(f>) - VfdW] ■ 

According to ( |2.13| ), h approaches a finite limit for A; ^ 0. Upon making the change 

of variables p — >■ —p — k the contribution proportional to h is seen to be cancelled by 
the first integral in (|2.14|) . We therefore conclude that 

limfl£'''H0,I^) = -iNjgHjm f ^h{(l),^P,r],^,g) A fiM'P,^) • (2.15) 



We have now dropped the subscript "f.T.", since in this limit only ( p.l4| ) contributes 
to the screening mass. To calculate this limit we proceed as before and set = (j) + e. 

One then verifies that for e — (or A; — > 0) the fermionic contribution to the screening 
mass squared is given by 

\ml(/3, a, m)] f^rm = Nrg^/J / { ^ + ^ ) . (2.16) 

L eiVM.P. )\}erm ^ (27r)3 Ug/^C-^+A) + 1]2 [e/3(</,-/i) ^_ 1]2 J ^ > 

In the continuum limit the electric screening mass is given by 

^li = lim ^iTT'lii-^ f^a, ma) . (2.17) 

For Wilson fermions, only momenta p in the immediate neighbourhood of p = 
contribute to the integral ( |2.1(j| ) for /3 — oo, = /3m = /3m fixed. But in 
this limit /30(p) —>■ (S^/p"^ + m?. Introducing in ( p.l6| ) the dimensioned momenta 
p = p/a as new integration variables, with a the lattice spacing, one then verifies, 
after performing the angular integration, and a partial integration that 

[mli]ferm = Nf— / dp [vfd{E, /i) + f]FD{E, /i)] , (2.18) 
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where E = ^/p^~l^m?, and 



r]FDiE,n) 



(2.19) 



are the Fermi-Dirac distribution functions for particles and antiparticles. 

We next consider the contribution to the electric screening mass arising from di- 
agrams (c-g) in fig. 1. They only involve sums over Matsubara frequencies of the 
bosonic type. Since in the continuum limit the only dimensioned scale is the tem- 
perature, their contribution to the screening mass will be of the form const x gT. 
For finite lattice spacing, however, the temperature dependence will be modified by 
lattice artefacts. In the following we first consider diagrams (c-e) which have an 
analog in the continuum. 

Hi) Contribution of diagram (c) 

Using the lattice Feynman rules given in Appendix A one finds after some algebra 
that this diagram contributes as follows to n44''^''(0, fc), defined in (12] 



3 ol 



mr(0,A;)(.) = V4 E 



TT (27r) 



where 



and 



a{k){z^ - If - h{q, k)z{z + lf 



a{k) = ^cos^^ 



b{q,k) 



(2.20a) 

(2.20b) 

(2.20c) 
(2.20d) 



Here p is generically defined by p = 2 sin The zeros of the denominator in (PT^UH) 



are located at 



where 



zi = ; Z2 = e ^ 
Z2 = ; z, = e-^ 



(j) = arcoshi7(g) , 
ip = arcoshi7(g + k) , 
H{P) = l + 2Esin^|. 



(2.21a) 



(2.21b) 
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The frequency sum can be calculated by making use of ( |B . 1| ) . After some straight 
forward algebra one finds that 



where 



and 



■h{(f),tp,a,b)fiBE{4>) 



^[h{4),i!,a,b) + h{^ip,4),a,b)] 



h{4>, ip, a, h) 



—a sinh^ (j) + |6[cosh + 1] 



sinh 0[cosh — cosh ip] 
1 



1 

is the lattice version of the Bose-Einstein distribution function. 
In obtaining this result we have made use of the fact that 



(2.22a) 

(2.22b) 
(2.22c) 



7 q-*-q~k ~ 
(p < > yj 



(2.23) 



while a{k) and b{q, k) are invariant under the transformation q —q — k. Note that 

the function ( |2.22b| ) is singular for k —>■ 0, since in this limit (p —>■ tp. The singularity 
is however integrable as can be seen by making use of ( 2.23 ) to write ( |2.22a ) in the 
form 



nS''^)(0,A;)(,) = 3<7 



:h{(j),'ijj,a,b) A fiBE{(t),i^) 



(2.24a) 



+ ^(?2|a(l!) + |_"^^^[/i(0,^,a,6) + /i(7/;,0,a,6)][l + 2f7BB(0)]| , 



where 



A t)be{<P, i^) = t1be{4>) - VbeW 



(2.24b) 



The limit k ^ can now be easily taken and one obtains the following contribution 
to the electric screening mass 



. ^ 2x 32/0 r 



3 coth 6 + 



2 sinh ( 



[i+2fjBEm 



15 r d^Q e' 
+ -ji-B 



13^ 



(27r)3 [eW - If 



(2.25) 
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iv) Contribution of diagram (d) 

This diagram involves the 4-gluon vertex, which consists of types of terms differ- 
ing in the colour structure: terms involving the structure constants /abC) and terms 
involving the completely symmetric colour couplings dABC- We denote the corre- 
sponding contributions to ti^fj^^\k) by [1144'^^ (A;)] [j] and [1144''^^ (A;)] [^j, respectively. 

Consider first [nit''^(0, After some algebra one finds that 



id) 



3 1 



(2.26a) 



where 



and 



[/] 



-c(fc)(^^ + 1) + dCk){z - 1)^ + P(g, k)z 
[z - Zi][z - Z2] 



1 + 2 ^ COS A; 



■3 1 



c{k) 
dil) 
Pill) 

Performing the frequency sum in ( p.26a| ) one obtains 

fif,'\o,k)^ 



1 ~2 

2 + lY.[i^Tk)] + iq~k)]] . 



+ 



id) =^g^{-c{k) + d(k) 
[f] 4 I 

''^ d^Q ~ 11 

^ (c-c/)coth0+ (rf- -P) — 



TT (27r) 



2 sinh I 



(2.26b) 



(2.26c) 



[l + 2r)B£(0)] . (2.27) 



Taking the limit A; — > one finds the following contribution to the screening mass 



d^q 



^"^^'^[/i = /r'+2 7-. (2vr)3 



Next consider the contribution \u!f^^\{]^k)(^d)\[d]- It is given by 



17coth(/) + 



sinh I 



[l + 27]BE{m ■ (2.28) 



Iift\Q,k) 



2 



[d] 



(2.29a) 
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where 



K{k) 



1 ^20 , K(k){z-lf-L{lk)z 
12 Vy + ^^^^J [z - -z,]\z - -z,] 



x2 



i 

and where d{A) is defined by 



(2.29b) 
(2.29c) 

(2.29d) 



(2.29e) 



E,F=1 

Performing the frequency sum one finds 



mr(o,^)(,) 



1 /20 

24 vy 



^7^7 ( =f + ^(^) 



(2^ 



ii' + -L 

Xcoth0 ^ 

sinh 6 



[1 + W(0)] 



Since i^(^) and L(g, ^) vanish for ^ — > 0, it does not contribute to the screening 
mass, i.e., 

[(meO(d)][d] = . (2.30) 

v) Contribution of diagram (e) 

The only other diagram possessing a continuum analog is the ghost loop shown 
in fig. le. Its contribution is given by 



where 



1 (z''-iy 



2UU[z-Zi] 
Performing the frequency sum one finds that 



(2.31a) 



(2.31b) 



n2''^^(o,^)(e) = y 



+ 4 



(2^ 



(2.32a) 
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where 



sinh 
cosh — cosh %l) 



(2.32b) 



This function is again singular for A; ^ 0. To compute the hmit we proceed as 
discussed earher and write ( p.32a| ) in the form 



+ 2 



. (27r; 



(27r)3 

5((0, V;) A TiBEi(i),i') 



where AfiBE{4>,'ip) has been defined in ( 2.24b|) . Taking the hmit k ^ one obtains 

3 



el){e) 



4^ i-' + 7-.(2vr)3 



-2/3 



[l + 2f/Bs(0)]coth. 



n (27r)3 _ i]2 



(2.33) 



Combining the resuhs ( |2.25| ), ( |2.28| ), ( |2.30| ) and ( p.33[ ), we therefore find that those 
diagrams possessing a continuum analog yield the following contribution to the elec- 
tric screening mass in the gluonic sector 



(m 



el){c+d+e) 



2J_3 , i 

^ [ 4^27-. (27r)3 



coth (h — 



1 



sinh I 



[1 + 2r}. 



BE 



2" i'^ d?q 



(27r)3 [e/3<^ - l]2 ■ 



(2.34) 



The computation of the remaining contributions arising from diagrams (f ) and (g), 
which are a consequence of the lattice regularization, is straight forward. One finds 
that they cancel the first term in ( |2.34| ). Hence the gluonic sector (G) contributes as 
follows to the screening mass 



n (27r)3 [eP<^ - If 



(2.35) 



In the continuum limit ( |2.17| ) the corresponding expression for the (dimensioned) 
screening mass squared is given by 



(m,jG = hmQg (3 



f (27r)3 [ei/3^W) _ 1]2 



poo 

dqq' 



[e^t - 1] 



2 ' 
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where q= \q\. After a partial integration this expression takes the form 

6 n_o r°° , X 



b f 
('^ei)G = -^9^T^ / dx- 

Making use of 

POO x'^~^ 

/ dx— — - = r(a)C(«) , a > 1 , 
JO — 1 

where T{a) is the Euler Gamma function, and C{<^) the Riemann Zeta- function, we 
recover the well known result: {m1i)G = g^T^ ■ 



3 Lattice Artefacts in the Screening Mass 

In this section we compare the one loop result for the electric screening mass on the 
lattice with the continuum. This will provide us with an estimate of the magnitude 
of the lattice artefacts to be expected in numerical simulations. The numerical data 
we present is for two mass-degenerate quarks. 

On a lattice the temperature can be varied by either keeping the lattice spacing 
fixed and varying the number of temporal lattice sites, or by varying the lattice 
spacing (or equivalently the coupling) keeping A^,- fixed. For fixed lattice spacing the 
dependence of the screening mass on the temperature, fermion mass and chemical 
potential is given by (see eq. ( |2.17|) ) 

1 1 

[mei{T,m, n,a)]iau = -mei{—^^ia,ma) (3.1) 
If the lattice expression is to approximate the continuum, then the lattice spacing 




Figure 2: Dependence of [mei\iatt/[mei\cont on — for /z = and different lattice spac- 
ings measured in units of m~^. Open squares (filled stars) correspond to lattices with 

Nr = S {Nr = 16). 

must be small compared to all physical length scales in the problem. Hence we must 
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have that a <^ j;, a <^ and a <^ K We therefore expect that for temperatures 



T ^ - the continuum is well approximated for ma <^ 1 and /ia ^ 1. This is shown 
in figs. 2 and 3 where we have plotted [mei\iatt/[mei]cont as a function of T/m at /i = 
and /i/m = 1.5 for various lattice spacings measured in units of m~^. For = and 
ma = a G [0, .3], the deviation of this ratio from unity is seen to be at most 1.75% for 
^ < Y^a' '^^^ '^^^ '-'^ ^^^^ inequality is the temperature associated with a lattice with 
16 sites in the temporal direction. For ^ ~ ^ the deviation is already 7.5%. The 
endpoint of the curves for a = 0.05, 0.1 and 0.3 correspond to the minimal number 
of temporal lattice sites, i.e. A'^^ = 2, and the open squares (filled stars) to Nr = 8 

(Nr = 16). 



[melllatt 
[mellcont 



a = 0.3 



a = 0.1 



a = 0.05 




T/m 



Figure 3: Dependence of [mei]iau/[fnei]cont on ^ for ^ = 1.5 and different lattice 
spacings measured in units of m~^. Open squares (filled stars) correspond to lattices 
with Nr = 8 {Nr = 16). The dashed-dotted (dashed) fines sfiow tfie temperature 
dependence for a fixed number of temporal lattice sites, Nr = 8 {Nr = 16). 



For fi/m = 1.5 we must also ensure that /xa -C 1. We therefore expect that the 
allowed range of lattice spacings for achieving an accuracy of 2% for T < is now 
restricted to a smaller interval. Fig. 3 shows that the continuum screening mass is 
well approximated for a < 0.1 in this temperature range. 

In lattice simulations one is interested in determining the electric screening mass 
above the deconfining phase transition. It is extracted from correlators of Polyakov 
loops [6-9], or from the long distance behaviour of the gluon propagator [|1^]. In 
these simulations the number of lattice sites Nr = /3 is fixed. The electric screening 
mass in physical units, divided by the temperature, is then given by 

[mei]iau/T = Nrm,i{Nr, {fi/T)N;\ {m/T)N;') , (3.2) 

while in the continuum limit {Nr — > oo) this ratio is just a function of m/T and /i/T. 
The above conditions for approximating the continuuum now read (a) A^,- 2> 1; (b) 
ma = {y)-^ ^ ^'^ = (t^W ^ therefore expect that for fixed A^,- the 

continuum is best approximated for high temperatures. For = 1.5 this is shown in 
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fig. 3, where the temperature dependence of [mei]iatt/["mei\cont for Nr = S and = 16 
is given by the dash-dotted and dashed curves respectively. The strong deviation of 
this ratio at low temperatures is due to the fermion loop contribution. This is evident 




Figure 4: Dependence of the pure gluonic contribution to [mei]iatt/[Tnei]cont on the 
number of lattice sites Nr. The solid line interpolates between different numbers of 
lattice sites, A^,-- 

from fig. 4, where we have plotted [mei]iatt / ['m'ei]cont for the pure SU(3) gauge theory 
for various values of Nt- This ratio only depends on the number of lattice sites, A',-. 
The solid line interpolates between different numbers of temporal lattice sites. The 
deviation from the continuum is seen to be small already for A^^- = 8. For A^^- = 8, 
and Nj. = 16, it is about 2%, and 0.4%, respectively. 




Figure 5: Dependence of [mei]iatt/['m'ei]cont on ^ and ^ for (a) A^,- = 8, and (b) 
A^^ = 16. 

Finally, in fig. 5 we have plotted the ratio [mei]iatt/[iTT'ei]cont as a function of ^ 
and ^ for A^,- = 8, 16. In the parameter range considered the above inequalities are 
well satisfied and the deviation of this ratio from unity is seen to be at most 1.7% 
for Nt- = 16, and 10% for A^^ = 8. The range of values m/T and fi/T for which the 
continuum is well approximated will of course increase with increasing A^.^- 
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4 Conclusions 



In this paper we have have computed the electric screening mass for Wilson fermions 
in the infinite volume limit for lattice QCD at finite temperature and chemical poten- 
tial in one loop order. The expression we obtained had a very transparent structure 
in which the artefacts arising from a finite lattice spacing were concentrated in two 
functions which in the naive continuum limit reduced to the on shell energies of a 
free quark and gluon. We then studied the dependence of the lattice screening mass 
on the temperature and chemical potential for fixed values of the lattice spacing. It 
was found that lattice artefacts give rise to an enhancement of the screening mass. 
For /i = and lattice spacings a < .3m~^ the deviation of [mei]iau/[^ei]cont from the 
continuum was found to be less than 1.75% for — < t^, where a = ma. For = 1.5 
a substantially smaller lattice spacing was required to approximate the continuum. 
Most of the deviation was found to be due to the fermion loop contribution. 

Since in numerical simulations the temperature dependence of the screening mass 
is extracted from a given lattice by varying the coupling (lattice spacing), we have 
also studied the dependence of the lattice screening mass on the temperature and 
chemical potential for fixed N^. It was found that for N^- = 16, and temperatures 
larger than the fermion mass and chemical potential, the continuum screening mass 
was approximated to within 1.75%. The corresponding deviation for A^,- = 8 was 
found to be at most 10% and to be due to the fermion loop contribution. In the 
pure SU(3) gauge theory the continuum was already approximated to 2% for only 8 
lattice sites in the temporal direction. 

Our analysis was carried out for infinite lattice volume. For finite spatial vol- 
ume the momentum spectrum becomes discrete and zero momentum modes must 
be treated separately in a perturbative expansion. In comparing the data for the 
electric screening mass obtained in Monte Carlo simulations with continuum pertur- 
bation theory, finite volume effects need to be included, while finite lattice spacing 
effects may, as our analysis suggests, be negligible for those couplings and lattice 
sizes at which the simulations have been performed. 
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A Feynman rules 

For vanishing temperature and chemical potential the lattice Feynman rules have 
been given in ||T^. The expression for the four gluon vertex stated in that reference 
is however not quite correct. For completeness sake we collect in this Appendix the 
Feynman rules for lattice QCD with Wilson fermions at finite temperature and chem- 
ical potential. In this case the fourth component of boson and fermion momenta are 
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replaced by k4 — uj'l — 2t:1/ (5 and p4 = ujj +1/1 — {2i + l)7r//3 + i// respectively. 
As usual the gluon fields are defined at the middle of the links. Following standard 
conventions, we did not include the combinatorial factor arising from the symmetriza- 
tion of the vertices. The dimensionless version of these rules, used in section 2, are 
obtained by setting a ~ 1. 



i) Fermion, gauge field and ghost propagators: 

7m sin(ap^)+m(p) 



b,/3 
where 



m 



(p) =m + ^^sin2(%^) 



is the momentum dependent Wilson mass and r the Wilson parameter which we have 
set to 1 in the calculation. 



^^n^W^ Sab^, [S,. - (1 - A) 
B,i/ k A,fj, 



B 



k 



Sab 



fe2 



ii) The vertices: 

>a(p) 




{27rrSl){-p + q + k) 



ig{T^)ab{{lix)ai3 cos(f (p + q)^) - ivSa/s sin(f (p + q)i,)} 



^^{2n)%,^{-p + q + k + f) 



-(55,, 

\ agH,,{T^, T^}ab{r5^^ cos(f (p + q)^) 
-i{l,,)a(} sin(f(p + g)^)} 
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/ 



igfABcfi^cos{^s^) 




PSo,^.^++^l+^l+^l) i27if6l,){-f+ s + k + q) 



12 



-(58,, 




WfABc{K\{k - r)^ cos(| g^) 
+8nv{q - k)^ cos(f r^) + 5x„{r - q)^ cos(f /ca)} 




where 



{2T:f5lS-rq-rr-rs) 



'-,(A),ABCD 



{k,q,r,s) 



rjfj^f ^■°(^' ^' ^' = I fABEfcDE 



- Sij,xSup{l2 cos{^a{k - r)^ cos{^a{q - s)^) - a^Kq^f^s^] 
+ 5^p5^A{12cos(^a(/c - s)^cos(^a(g - r)^) - a^kuquf^^s^} 
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and 



+ 5i,x6^p2o?{s - r)^qf, cos(^aA;j.) 

- S^J^p2a'^{q - k)^rpCOs{^asx) 
+ S^Ji^x'^a^iq - k)pSx cos( ^arp) 

- S^,J|,xS^,pa^ - k)^{s - r)Jj 

+ the cyclic permutations: 
/ {B,q,iy) (C,r,A) {D,s,p) \ 
[ (C,r,A) iD,s,p) iB,q,u) ) 

• {B,q,iy) {C,r,X) {D,s,p) \ \ 
iD,s,p) iB,q,u) (C,r,A) Jf 

+ 1-{SabScd + SacSdb + SadSbc) 

+ y^XdABEdcDE + dACEdoBE + dADEdBCE) \ 
E ) 

^IJiv5^x5,jip XI KQafaSa + ^ ^vhpkxQxf ,j.S p, 

+ ^iiX^upkuQ^jTv^p, + 5 iip^uxk'uqp'f' — (^uX^^upkuq^fpSp 
-5^ix5p.pK%fyS^ - Sp^ySp^pkxqxfpSx - Sp,^6pxkpqpfpSp,}} . 

B Frequency summation formulae 

In this Appendix we derive summation formulae over Matsubara frequencies useful 
for performing calculations on a lattice, where the frequencies are restricted to the 
Brillouin zone. We will first consider the bosonic case and afterwards turn to the 
fermionic one. 

B.l Bosonic frequency sums 

We will prove the following summation formula: 

Let g{z) be a meromorphic function of a complex variable z that is bounded 
for 1 2; I — > 00 and has no singularities on the circle \z\ — 1. Then 

1 



p 



n=-$/2 



Z, 



where /3 is a positive even number, o)^ = 21:71/(3 and Res^. {^g{z)^ denotes 
the residue of -g{z) at the pole Zi. 
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Consider the function 



h{u;) 



ip 



exp{i(3od) — 1 

which has simple poles at = lj^, n & Z, with residue 1. Because of the above 
conditions on g{z) there exists an e > such that g{e^^) has no singularities for 
Im(c2;) G [— e,e]. Hence 



/3/2-1 



(B.2) 



where C is the integration contour depicted in fig. B.l. Introducing the variable 



C 



^ I- 



-TT 



TT 



Figure B.l: Integration contour C in eq. ( p.2| ) 

z = exp(za)) we obtain 



,3/2-1 



n=-/3/2 



27rz Jc, 



(B.3) 



where Ci and C2 are the contours shown in fig. p.2| with the integration carried out 
in the indicated sense. 











— 


Ci \ \ C2 


\-l \ 


/ 1/ 











Figure B.2: Integration contours Ci and C2 in eq. 
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Since limi^i^oo z~^g{z) = we can distort the outer integration contour to infinity 
taking into account the singularities of g{z) for \z\ > 1. In the case where g{z) is a 
meromorphic function ( B.l ) follows immediately. 



B.2 Fermionic frequency sums 

In the fermionic case we will proceed in an analogous way to prove the following 
summation formula: 

Let g{z) be a meromorphic function of a complex variable z that is bounded 
for \z\ — > oo and has no singularities on the circle \z\ = exp(— /x). Then 



I /3/2-1 
1^ m=-/3/2 



E 



Res,, {Igjz)) 
e^Hf + 1 



(B.4) 



where Cj^ = (2m + l)vr//3. 
Consider the function 

exp(i/?(u; — i/t)) + 1 

which has simple poles at cD = ci)~ + ifi, m E Z, with residue 1. Because of the 
conditions on g{z) there exists an e > such that g{e^^) has no singularities for 
Im(ci)) G [/i — e, /i + e]. Hence 



/3/2-1 



1 



^ g^S'^-^m) = - — jjoj g 



m=-/3/2 
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exp(i/3(ci) — i/i)) + 1 



(B.5) 



where the integration contour C is depicted in fig. [B.3|. Introducing again the variable 



C 



• • • • * 



-TT 



TT 



Figure B.3: Integration contour C in eq. ( p.5| ) 

z = exp(za)) we obtain 



I /3/2-1 
m=-/3/2 



2ni^Jc, z e'^f'zP + l 
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where the integration contours Ci and C2 are as depicted in fig. |B.2| except that 



the radii are changed to e~^~^ and e"^"'"'^ respectively. Hence, for the same reason 
as in the bosonic case, we can distort the contour C2 to infinity and are led, for a 
meromorphic function g{z), to the summation formula ( p.4|) . 
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